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ABSTRACT 

An optimal-observable analysis of the angular and energy distributions of the 
leptons and bottom quarks in the process e+e~ — > — i^/ b ■ ■ ■ has been per- 
formed in order to measure the most general top-quark couplings to gauge bosons 
at polarized linear colliders. The optimal beam polarization for determination of 
each coupling has been found. A very sensitive test of CP violation in ti production 
and decay has been proposed. 
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1. Introduction 

In spite of tlie fact that the top quark has been discovered already several years ago 

its interactions are still very weakly constrained. It remains an open question 
if top-quark couplings obey the Standard Model (SM) scheme of the electroweak 
forces or there exists a contribution from physics beyond the SM. We could interpret 
the great success of the 1-loop precision tests of the SM as a strong indication that 
the third generation also obeys the SM scheme. However, an independent and direct 
measurement of the top-quark couplings is definitely necessary before drawing any 
definite conclusion concerning non-standard physics. 

Over the past several years there was a substantial effort devoted to a possibility 
of determining top-quark couplings through measurements performed at the open 
top regionQ of future e~^e~ linear colliders The existing studies focused 

mainly on tests of CP violation in top-quark interactions. In this article we will 
construct some new tools which could help to measure both CP violating and CP 
conserving top-quark couplings at linear colliders and therefore reveal the structure 
of fundamental interactions beyond the SM. 

The top quark decays immediately after being produced and its huge mass 
rrit ~ 174 GeV leads to a decay width Pf much larger than vIqcd- Therefore the 
decay process is not influenced by any fragmentation effects and decay products 
will provide useful information on top-quark properties. Here we will consider 
distributions of either in the inclusive process e^e~ ^ ti ^ ■■■ or bottom 
quarks from e+e^ tt ^ b ■ ■ ■■ It turns out that the analysis of the leptonic and 
6-quark final states is similar and could be presented simultaneously. Although ti 
are also produced via WW fusion , we do not consider here this mechanism since 
a{e^e~ — ttuu) is expected to be much smaller than a{e^e~ — > tt) for the energy 
of our interest (v^ < 2 TeV) g. 

This paper is organized as follows. First in sec. 2 we describe the basic frame- 
work of our analysis, and then show the angular and energy distributions of the 

"^Recently an interesting and complementary analysis by Jezabek, Nagano and Sumino has 
been published |^ where the authors discussed possibility of determining CP- violating production 
form factors at the ti threshold region. 



lepton and 6-quark in sec. 3. In sec. 4, after briefly reviewing the optimal-observable 
procedure |T^ , we estimate to what precision all the non-standard parameters can 
be measured or constrained adjusting the initial beam polarizations. Finally, we 
summarize our results in sec. 5. In the appendix we collect several functions used 
in the main text for completeness, though some of them could also be found in our 
previous papers @, . 

2. Framework and Formalism 

We parameterize ti couplings to the photon and the Z boson in the following way 
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denote the SM contributions to the vertices. Among the above non-SM form fac- 
tors, 6A^, 6B^, 6Cv describe CP-conserving while (5D„ parameterizes CP- violating 
interactions. Similarly, we adopt the following parameterization of the Wtb vertex 
suitable for the t and i decays: 
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where Pl/r = (1 =F75)/2, Vtt is the (tb) element of the Kobayashi-Maskawa matrix 
and k is the momentum of W. In the SM //" = //' = ! and all the other form 
factors vanish. On the other hand, it is assumed here that interactions of leptons 
with gauge bosons are properly described by the SM. Throughout the calculations 
all fermions except the top quark are considered as massless. We also neglect terms 
quadratic in the non-standard form factors. 



Using the technique developed by Kawasaki, Shirafuji and Tsai [11] one can 
derive the following formula for the inclusive distributions of the top-quark decay 



product / in the process e+e ^ —*>/ + ■■ ■ [||: 
(e+e / + ■■■) =4 / dr2t—(n,0) 



where Ft is the total top-quark decay width and rf^-T/ is the differential decay rate for 
the process considered. da{n, 0) / df2t is obtained from the angular distribution of ti 
with spins s+ and s_ in e"^e~ tt, da{s^, sJ)/dQt) by the following replacement: 
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^ 0, (2.4) 
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where the matrix element for t(s+) —>/-!-■•■ was expressed as But{pt, s+), /!+ = 
+ rrit, d$ is the relevant final-state phase-space element and J2 denotes the 
appropriate spin summation. 

3. Angular/Energy Distributions 

In this section we present d^a/dXfdcos6f for the top-quark decay product /(= 
i^/ b ), where Xf denotes the normalized energy of / defined in terms of its energy 



Ef and the top-quark velocity /3(= Jl — Ami/ s) as 



2Ef l-P 



mt V 1 + 

and Of is the angle between the e~ beam direction and the / momentum, all in the 
e~^e~ CM frame. 

Direct calculations performed in presence of the general decay vertex ( p.2| ) lead 
to the following result for the n^^ vector defined in eq.( |2.4| ): 

% = an g^^ — Pf (3.1) 

V mi ) ptPf 

where for a given final state /, is a calculable depolarization factor 

1 for / = £+ 

H 2^ 8VF(1-.) , (3-2) 

2r + lL (2r-l)(2r + l) ^''^ ' i 



with r = [Mw/fTttY- Similarly we have = —a^ with replacement — > /I". It 
should be emphasized here that the above result means that there are no corrections 
to the "polarization vector" for the semileptonic top-quark decay. On the other 
hand, one can see that the corrections to a* could be substantial as the kinematical 
suppression factor in the leading term 2r — 1(= —0.56) could be canceled by the 
appropriate contribution from the non-standard form factor f^. 

Applying the strategy described above and adopting the general formula for the 
ti distribution da{s^, sJ)/dfIt from refs.|^, 12|, one obtains the following result for 



the double distribution of the angle and the rescaled energy of / for longitudinally 
polarized e^e~ beams: 

where is the fine structure constant and Bf denotes the appropriate branching 
fraction. The energy dependence is specified by the functions 6>/**'(xy), explicit 
forms of which for unpolarized beams were shown in ref. [|T^.0 They are parame- 
terized both by the production and the decay form factors. 

The angular distribution for / could be easy obtained from eg. ( p. 3D by the 
integration over the energy of /: 



d cos 0, 



r dx, /f = '-^^B, + Q[^'^ cos9, + Qi^'^ cos' 9, 
Jx. dXfd COS Of 2s ^ 

(3.4) 

where = /^^ (ix©/'*' are shown by eq.( |A.l| ) in the appendix and x± define 

kinematical energy range of x: 

r{l- p)/{l + p) <x,<l and (1 - r)(l - /5)/(l + /?) < x,, < 1 - r. (3.5) 



The decay vertex is entering the double distribution, eq.( p.3| ), through i) the func- 
tions F^{Xf), G^{Xf) and -fff 2(2^/) defined in the appendix, and ii) the depolarization 
factor . All the non-SM parts of F^, and if/ 2 disappear upon integration over 

"^Thc functions 0l''*\xf) for polarized beams could be easily obtained from formulas for unpo- 
larized beams replacing Dv.a,vAi Ev,a,va, -Fi~4, Gi^4 defined by eg. ( A. 18 ) with Dy^j^ yj^, Ey\ y^, 



Fj\^4, G\^^ as in eq.(A.17) in the appendix 



the energy Xf both for C'^ and 6, as it could be seen from the exphcit forms for i?/**'. 
Since = 1 for the leptonic distribution, we observe that the total dependence 
of the lepton distribution on non-standard structure of the top-quark decay vertex 
drops out through the integration over the energy [131."^ However, one can expect 
substantial modifications for the bottom-quark distribution since corrections to 
could be large. 

The fact that the angular leptonic distribution is insensitive to corrections to 
the V — A structure of the decay vertex allows for much more clear tests of the 
production vertices through measurements of the distribution, since that way we 
can avoid a contamination from a non-standard structure of the decay vertex. As 
an application of the angular distribution let us consider the following CP- violating 
forward-backward charge asymmetry:^ 
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where Pe- and Pe+ are the polarizations of e and e beams, (io"+/~'^*) is referring to 
/ and / distributions respectively, and Cm expresses the experimental polar-angle 
cut. As 6f ^ 71 — 9 J under CP, this asymmetry is a true measure of CP violation. 
Since da^^*^ / dcosOj is obtained from da^^*y d cos 9f by reversing the sign of cos^^ 
and Fi*} terms and replacing with —a-'' in i7o'i,2 5 the asymmetry is explicitly 
given by the following formula 
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with (in the leading order) 
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"■^The same conclusion has also been reached through a different approach using the hclicity 
formalism in ref . JT^ . 

Which is an integrated version of the asymmetry we have considered in rcf . p3[ . 
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where all the coefficients are specified in the appendix, the subscript (0) indicates 
the SM contribution and we expressed as + a{Re{f^) with 
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As one could have anticipated, the asymmetry for f = £ is sensitive to CP 
violation originating exclusively from the production mechanism: It depends only 
on FiJ that contains contributions from the CP- violating form factors 6D^ and 6Dz 
while the contributing decay-vertex part consists of SM CP-conserving couplings 
only. For bottom quarks the effect of the modification of the decay vertex is 
contained in the corrections to b and b depolarization factors, a^ + a^ = a\lie{f^ — 
/^), with SM CP-conserving contributions from the production process.0 

It will be instructive to give the following remark here: The asymmetry is 
defined for various initial beam polarizations Pe±. For Pg- ^ Pe+, the initial state 
seems not to be CP invariant and therefore one might expect contributions to the 
asymmetry originating from the CP-conserving part of the top-quark couplings. 
However, as it is seen from eg. ( p. 81) , this is not the case. It turns out that even for 
Pg- 7^ Pe+ the asymmetry is still proportional only to the CP-violating couplings 



"^One can show that f^'^ — ifi'^ and /^'^ — ±/^'^ where upper (lower) signs are those for 
CP-conserving (-violating) contributions p5[ |. Therefore, when only linear terms in non-standard 
form factors are kept, any CP- violating observable defined for the top-quark decay must be 
proportional to f^'^ — fi'^ or /^'^ — /^''^. 



embedded in -Fi,2,3,4- The explanation is the following: Whatever the polarizations 
of the initial beams are, the electron (positron) beam consists of e(±l)(e(±l)) 
where ±1 indicates the helicity, and only e(±l) and e(=i=l) can interact non-trivially 
in the limit of rrie = since they couple to vector bosons. Therefore the interacting 
initial states are always CP invariant. 



Now, since we have observed in ref.|13] that the differential version of the asym- 
metry discussed here could be substantial for higher collider energy, in order to 
illustrate the potential power of the asymmetry we present in tabs.|l] and ^ (as a 
function of y/s) the expected statistical significance (Nsd) for the asymmetry: 



(1) Pe- = Pe+ = 



x/i (GeV) 


500 


700 


1000 


1500 




-1.2 ■ 10-2 


-2.6 ■ 10-2 


-4.0 ■ 10-2 


-5.4-10-2 




2.42 


3.87 


4.32 


3.94 



(2) Pe- = Pe+ = +0.8 



(GeV) 


500 


700 


1000 


1500 




-1.4-10-2 


-2.6 - 10-2 


-3.6-10-2 


-4.2 - 10-2 


Nso 


2.80 


4.09 


4.04 


3.25 



(3) Pe- = Pe+ = -0.8 



(GeV) 


500 


700 


1000 


1500 


A" 


-1.2 - 10-2 


-2.6 - 10-2 


-4.1 - 10-2 


-5.7-10-2 


Nso 


3.53 


5.72 


6.57 


6.20 



Table 1: The CP- violating asymmetry A^^p and the expected statistical significance 
Nso for Re{5D^^z) = +0.05, and beam polarizations Pe- = Pe+ = (1) 0, (2) +0.8 
and (3) —0.8 as an example. 



(1) Pe- = Pe+ = 



v^(GeV) 


500 


700 


1000 


1500 




+1.2 ■ 10-2 


+1.7-10-2 


+2.2 ■ 10-2 


+2.6 ■ 10-2 


Nsn 


5.10 


5.50 


5.03 


4.03 



(2) Pe- = Pe+ = +0.8 



v^(GeV) 


500 


700 


1000 


1500 


A'' 


-9.4 ■ 10-3 


-4.6 ■ 10-3 


+ 1.4-10-3 


+7.8 ■ 10-3 


Nsn 


4.04 


1.52 


0.33 


1.27 



(3) Pe- = Pe+ = -0.8 



v^(GeV) 


500 


700 


1000 


1500 


A*" 


+2.6 ■ 10-2 


+3.0 ■ 10-2 


+3.3 ■ 10-2 


+3.5 ■ 10-2 


Nsn 


16.0 


14.3 


11.2 


8.02 



Table 2: The CP-violating asymmetry A''^p and the expected statistical significance 
Nsn for Re{6D^ z) = — /|') = +0.05, and beam polarizations Pg- = Pe+ = 

(1) 0, (2) +0.8 and (3) -0.8 as an example. 



where Lg// = eL is an effective integrated luminosity for the tagging efficiency 

e. Hereafter we adopt the integrated luminosity L = 500 fb-^ and the efficiency 

e = 60% both for lepton and 6-quark detection.[^ In addition, to fit the typical 

detector shape we impose a polar-angle cut |cos6'^| < 0.9, i.e. = 0.9 in 

eq. (|3.8| ), both for leptons and bottom quarks. On the other hand, we will not 

impose any cut on the lepton/ 6-quark energy since their kinematical lower bounds 

^min 5 g^^^ ^rnin ^ 27.5 GcV (for ^ = 500 GeV) are large enough to be 

'^That low efficiency is supposed to take into account cuts necessary to suppress the back- 
ground. If the 6-tagging is apphed then, as shown in the second paper of ref. [^|, the irreducible 
background to top events due +26 + 2 j is negligible, provided that a vertex tagging effi- 
ciency et ^ .5 can be achieved. Therefore for the 6-tagging case the efficiency we have employed 
is definitely conservative. Since Nsd scales as VeL it would be easy to estimate the statistical 
significance for any given luminosity and efficiency. 



detected. Perfect angular resolution will be assumed both for lepton and 6-quark 
final states. Also ideal leptonic-energy resolution will be used. 

As it is seen from the tables, the asymmetry A^^^p turned out to be a very 
sensitive CP-violating observable; even for unpolarized beams and CP-violating 
couplings of the order of 0.05 one can expect 2.4cr ~ 5.5a effect both for lepton 
and 6-quark asymmetries once L = 500 fb~^ is achieved. 

The CP-violating form factors discussed here could be also generated within 
the SM. However, it is easy to notice that the first non-zero contribution to 6D^^z 
would require at least two loops. For the top-quark decay process CP violation 
could appear at the one-loop level, however it is strongly suppressed by the double 



GIM mechanism [|18]. Therefore we can conclude that an experimental detection 
of CP- violating form factors considered here would be a clear indication for physics 
beyond the SM. In particular, non-vanishing in the lepton distribution will 
strongly indicate some new-physics in tt'j/Z couplings. 

4. Optimal-Observable Analysis 
4.1. Optimal observables 

Let us briefly recall the main points of the optimal-observable (00) technique [T^ . 
Suppose we have a distribution 

^(=r(0)) = 5:Q/.(0) (4.1) 

where /j(0) are known functions of the location in final-state phase space and q's 
are model-dependent coefficients. The goal would be to determine Cj's. It can be 
done by using appropriate weighting functions Wj(0) such that / d(f)Wi{(f))IJ{(l)) = Ci. 
Generally, different choices for Wi{(f)) are possible, but there is a unique choice so 
that the resultant statistical error is minimized. Such functions are given by 



w. 



(0) = 5:X,,/,(0)/r(0), (4.2) 



where Xij is the inverse matrix of Mij which is defined as 



M,.id^mm. (4.3) 



The statistical uncertainty of Q-determination through da/dcj) measurement be- 
comes 



Aci = ^X,iaT/N , (4.4) 

where ctt = / d(f){da / dcj)) and N is the total number of events. 

It is clear from the definition of the matrix Mjj, eq. (|4.3| ), that Mij has no 
inverse if the functions /«(</>) are linearly dependent, and then we cannot perform 
any meaningful analysis. One can see it more intuitively as follows: if /i(0) = /j (0) 
the splitting between Cj and Cj would be totally arbitrary and only Cj + Cj could be 
determined. 
4.2. For application 

In order to apply the 00 procedure to the processes under consideration, we have to 
reexpress the distributions in the form shown in eq.( [4.1|) . The angular distribution, 
eq. (|3.4| ), has already an appropriate form for this purpose, where /j(0) = cos* 6'^ 
{i = 0, 1, 2) and i7/'*' are the coefficients to be determined. On the other hand, the 
double angular and energy distribution eq. (|3.3|) must be modified. We reexpress 
the distribution in the following way, keeping only the SM contribution and terms 
linear in the non-standard form factors: 



dxfd cos 9f 2 s 

where 



-B,S}*\x„e,), (4.5) 



+ [MSA.)J'!,':\xf,e,) + Re{SB,)H':\x„e,: 

v='y,Z 

+ Rei5C.)J'lj:\x,, e,) + Re(5D,)^i,W(x„ 6,) 



As it is seen from the above formula, the coefficients q of eq. (|4.1|) are just the 
anomalous form factors to be determined. The SM contribution reads: 

5f *)(a;,, e,) = Or\x,) + cos^, Or\x,) + cos' 9, 0^r'\xt) (4.6) 
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+ 2alRe{D\ 

Explicit forms of the functions ^j^^*^ c ^} and J^l^^ are shown in the appendix 
together with the functions f-'ix), g^ix) and h{ 2ix). 

There are ten functions entering eg. ([4.5]): Sf''*\ c d}{7 z} ■ 

explained earlier, one cannot determine their coefficients separately if they are 
not independent. As could be found from the appendix, for the double lepton 
distribution, the first nine functions are linear combinations of 

fix), fix) cose, f{x)cos^e, 



g\x), g\x)cos9, g\x)cos^9, 



/i'i2(a;)(l-3cos2 



/i^(x) cos 6', 



(4.10) 



while the last one, ■, is a combination of 5{f^,g^,h\2\ix) and cos"^ (n = 
0, 1,2). Since there are ten coefficients to be measured^ it looks always possible 
to determine all of them. However, it turns out not to be the case in some special 
cases. Indeed the possibility for the determination of all the ten form factors 
depends crucially on the chosen beam polarization. 

This can be understood considering the invariant amplitude for ee ti, which 
could be expressed in terms of eight independent parameters as 

Miee -^tt) = Cvv [Veli^Ue ■ Ua^Vt] + CvA [Velf^Ue " ^tTsT^^t ] 



'^Counting the SM coefficient in front of sf''*'' which is normahzed to 1 



+ CaV [VelblpiUe ' U^^Vt ] + C AA [VelblfjiUe " Ut'^bl^Vt ] 
+ CvS [ VeflUe " UfVt ] + Cyp [ Ve§Ue ■ Ua^Vt ] 

+ Cas [velbUUe " UtVt ] + Cap [ t'eTs^We ■ ua^vt]- 

However, if e or e is perfectly polarized, contributions from [ve'~i^u^ and [ve'~i^'~i^u^ 
are identical. For example, when e has fully left-handed polarization, u^. is replaced 
with UeL = (1 — 75)^e/2 and in that case they are changed as 

Therefore, the invariant amplitude becomes 
M{ee tt) 

= {Cvv + Cav) [Velf^UeL ' Mt7^^i] + {CvA + C aa) [Vel^iUeL ' Ut'^bl^Vt] 
+ {Cvs + Cas) [VeflUeL ■ UtVt] + {CvP + Cap) [VeUUel ■ Mt75^^t], 

and one ends with just four independent functions and therefore only four coeffi- 
cients could be determined. More details could be found in the appendix below 
eq.( |A.16D . Of course, such singular configurations of polarization are not considered 



in our analyses. 

As for 6-quark distributions 5f^{x) = 5g^{x) = 6h\{x) = dh^i^x) = 0, instead of 
ten functions we have in that case only nine of them given by the 6-quark 

version eq.( [4.10|) . Therefore, at most nine couplings could be determined. Since 
6-quark energy resolution is expected to be relatively poor, we will not apply 00 
procedure to the 6-quark double distribution. 
4.3. Numerical analysis 

Below, we will adjust beam polarizations to perform the best measurement of 
the form factors. In order to gain some intuition we show in figs. ^ and the 
functions ^[a'b c D}{-y z}^ -^m P^^^ Sf'*^ for unpolarized beams (fig.|I|) and for the 
beam polarization Pg- = Pe+ = +0.5 (fig.0). The figures illustrates how much the 
polarization could modify the functions and therefore infiuence the possibility for 
the determination of the form factors. 




Figure 1: The shape of the coefficient functions ^{ABCDyi-yZ}^ ^2r\ 
unpolarized beams 



Lepton angular distribution 

Since we have only three independent functions {1, cos 9, cos^ 9}, M and its inverse 
X are (3, 3) matrices. We have considered the following polarization set-ups: Pe- — 
Pf.+ — 0, ±0.5 and ±1. Since 1 > | cos 9\ > cos^ 9 we observe that Xu < X22 < X^s, 

(*) A (*) 

therefore the statistical uncertainty for measurement, AHq, is always the 
smallest one. 

Once we assume the detection efficiency e and the integrated luminosity L, we 
can compute the statistical significance of measuring the non- SM part of 

For the efficiency and luminosity specified earlier we obtain 



• Pe- = Pe+ = +0.5 



• Pe- = Pe+ = +1 



• Pe- = Pe+ = -0.5 



Mn 


= 2.06, 


M22 = 


- 0.55, 


M33 


= 0.27 




= 1.09, 


X22 — 


1.83, 




= 8.40 




= 16.1, 




= 3.5, 




= 1.9 


Mn 


= 3.06, 


M22 = 


: 0.95, 


M33 


= 0.49 




= 0.82, 


X22 = 


1.66, 


-^33 = 


= 6.25 




= 7.5, 




2.7, 




= 1.3 


Mn 


= 3.20, 


M22 = 


: 1.25, 


M33 


= 0.72 


Xn 


= 1.00, 


X22 = 


2.39, 


-^33 = 


= 7.00 




= 5.2, 




3.0, 


7V(2) . 


= 1.3 


Mn 


= 1.27, 


M22 = 


: 0.35, 


M33 


= 0.17 


Xn 


= 1.81, 


X22 — 


2.91, 


-^^33 = 


= 13.4 




= 26.2, 


Ar(i) = 


= 4.9, 




= 3.0 



• Pe- = Pe+ = -1 

Mil = 0.76, M22 = 0.21, M33 = 0.11 

Xn = 3.06, X22 = 4.90, X33 = 22.3 

7VW=35.5, 7Vi^)=6.4, TV^') = 4.0, (4.11) 

where we put all the non-SM parameters Ke{5{A, B, C, Dj^^z) and Re(/^) to be 
+0.05 as an example. 

As one can see, the precision is better for negative beam polarization, partly 
because of larger number of events. However we cannot conclude that using 
negatively-polarized beams is always more effective for new-physics search, since 
Nso strongly depends on the non-SM parameters used in the computations. In 
fact, positively-polarized beams give smaller and this is independent of the 
choice of non-SM parameters. Therefore polarization of the initial beams should 
be carefully adjusted for each tested model in actual experimental analysis. 
6-qu£Lrk anguleir distribution 

We can compute M, X and Nsd in the same way as for the lepton distribution: 

• P,- = Pe+ = 



• = = +0.5 



• Pe- = Pe+ = +1 



Mil 


= 2.23, 


M22 = 


- 0.63, 


M33 = 


= 0.31 


Xn 


= 1.04, 


X22 — 


1.85, 




7.98 




= 37.3, 


A^(^) -- 


= 17.8, 




= 1.9 


Mil 


= 2.38, 


M22 = 


: 0.65, 


M33 = 


= 0.32 


Xn 


= 0.96, 


X22 = 


1.61, 


-^33 = 


7.29 




= 15.5, 


-- 


= 7.6, 


7V(2) . 


= 1.8 


Mil 


= 1.63, 


M22 = 


: 0.45, 


M33 = 


= 0.22 


Xn 


= 1.39, 


X22 — 


2.29, 


-^^33 = 


10.5 




= 9.7, 


Ar(i) = 


4.8, 




2.2 



• Pe- = n+ = -1 



1\ /f 

Mil 


1 /I 

= i.4o, 


M22 


n /1 
= U.4Z, 


/I /f 


= U.zi 


Xn 


= 1.63, 


X22 = 


= 3.01, 


-^33 = 


12.5 




= 62.6, 


^ ^SD 


= 29.3, 




= 2.4 


Mil 


= 0.87, 


M22 -- 


= 0.25, 


M33 = 


= 0.13 




= 2.74, 


X22 = 


= 5.10, 


-^33 = 


21.0 




= 85.1, 


iV(i) 

^SD 


= 39.7, 


]V(2) 


= 3.1, 



for Re((5{y4, B, C, D}^^z) = Re(/|^) = +0.05. Negatively-polarized beams give 
better precision again, but the same remark as to the lepton angular distribution 
should be kept in mind also here. 

The above results prove that the optimal observables utilizing the angular dis- 
tributions should be very efficient seeking for the non-SM parts of i?/'*'. However, 
since they are combinations of the form factors, we can only constrain them. Of 
course, it would be exciting if we found any signal of non-standard physics, however 
our final goal is to determine each form factor separately. That is why we proceed 
to the next analysis using the double angular and energy distributions. 
Lepton angular and energy distribution 

Because of high precision of direction and energy determination of leptons we 
adopted the double energy and angular distributions, eq.( [4.5|) , also for 00 analysis. 
As discussed earlier, in principle all nine form factors could be determined with the 
expected statistical uncertainties Aci for q = Re{6{A, B,C, D}^^z) and Re(/^). 
The beam polarizations Pe- and Pe+ were adjusted to minimize the statistical error 
for determination of each form factor. We found that positive polarizations lead 
to a smaller Aci for eight form factors in the production vertices. Unfortunately, 
however, the optimal polarizations for each form-factor measurement is different. 
Below we present the smallest statistical uncertainties and the corresponding beam 



polarizations for each parameter: 



A 



A 



A 



A 



A 



A 



A 



A 



Re(M^) ] 

/ M^iO.lS, 

V 5Cz: 4.69, 

Re(Mz) ] 

/ M^: 0.23, 

V SCz- 1.76, 

Re(5B^) ] 

/ SA^: 0.43, 

V &Cz- 1.17, 

Re(dB^) ] 

/ SA^: 0.25, 

V SCz- 1.56, 

Re(5C^) ] 

/ SA-y-. 0.82, 

V 5Cz: 1.11, 

Re(5Cz) ] 



= 0.16 

<5B^: 0.25, 
5D^: 27.2, 

= 0.07 

0.11, 
SD-,: 7.09, 

= 0.09 

Mz; 0.11, 
5D^: 0.95, 
= 0.27 

6Az. 0.07, 



for P,- = 0.7 and P,- 



5Dz:53.2, /f:0.02 

for Pe- = 0.5 and Pe+ 

<5Bz: 0.27, SC^: 0.70 
SDz:20S, /f:0.02 



0.7, 



0.4, 



for P, 



0.2 and P,+ = 0.2, 



5Bz: 0.36, 5C^:0.21 



for P.- =0.4 and P,- 



0.11 



(5Dz:18.5, /I*: 

for P,- 



0.02 



0.1 and P, 



5Az;0.22, 5B.^;0.10, 5Bz:0S5 
5D-,: 1.76, 
= 1.11 



5Dz:14.6, /I*: 0.03 



for P, 



0.1 and 



/ dAj-.D.tiz, Q/lz;u.zz, 0£>-^;u.iu, 0£>z'- 
V (5C.^:0.11, 5D^:1.76, <5Dz: 14.6, 0.03 

Re((5D^) ] = 0.08 for = 0.2 and 

/ M^:0.52, Mz;0.13, 5B^;0.09, 5Bz:0A2 
\ 5C^f:0.15, 5Cz: 1.13, <5Dz: 14.4, f^-.O.OS 

Re{SDz) ] = UA for Pe- = 0.2 and Pe 

/ M.y:0.52, Mz; 0.13, 5B^:0.09, 5Bz:0A2 
\ (5C^:0.15, 5Cz: 1.13, 50.^:0.08, 0.03 



0.4, 



0.0, 



0.0, 



0.1, 



0.1, 



(4.13) 



where we also showed the expected precision of the other parameter measurements 
for the same beam polarizations. For instance, we can expect Z\[Re((M^) ] = 0.16 
for Pg- = Pe+ = 0.7 while the expected precision of Re((M^), Re((5B^), ■ ■ ■ for the 
same polarizations are 0.13, 0.25, ■ ■ ■, respectively. This result is independent of 
the choice of the non-SM parameters in contrast to the preceding results. 

As it is seen the precision of S{C, D}z measurement would be very poor even 
for the optimal polarization. This is mainly a consequence of the size of ^j^'/)}^, 
which is illustrated in figs.|l] and 0: These two functions are very small in a large 
area. More quantitatively, the size of the elements of M matrix, Mij, is 0(1) 
for i,j ^ 7,9, while the size of M„{= M-ji) and Mi9(= Mg^) is at most 0(10^2)^ 



In addition, determination of 6D^ would be practically difficult, as well, since its 
error varies rapidly with the polarization. For example, Zl[Re(5D^) ] becomes 0.86 
for Pg- = 0.1/Pe+ = 0.1 and 0.99 for P^- = 0.3/Pe+ = 0.1. The source of that 
sensitivity is hidden in the neutral-current structure with sin^ 9w — 0.23. Indeed, 
the optimal polarization becomes Pg- — 0-1 instead of 0.2 (Z\[Re(5D^) ] = 0.09) for 
sin^^^vi^ = 0.25. On the other hand, a good determination (almost independently 
of the polarization) could be expected for Indeed, the best precision is 

A[ Re(/2^) ] = 0.01 for P^- = -0.8 and = -0.8 (4.14) 

whereas even for the unpolarized beams we obtain A[ Re{f^) ] = 0.03. 

At the time a linear collider will be operating, data from Tevatron Run II and 

LHC will also provide independent constraints on top-quark couplings. Below we 
provide an example of a combined analysis assuming SA^, 6By and are known 
and 00 are used to determine 5C„ and 5D„ only (here we put SAy — SBy — — 
for simplicity). The results are as follows: 

A[Re{6C^)] = 0.04 
A[Re{6Cz)] =0.23 
A[Re{5D^)] = 0.03 
A[Re{5Dz) ] = 2.97 

The error for SDz became much smaller but still too large for practical use. How- 
ever, as we have seen in sec. 3, the CP-sensitive asymmetry A^p would provide much 
stronger constraints on SD^^z- 

5. Summary and Conclusions 

(-) 

We have presented here the angular and energy distributions for / in the process 

_ (-) 

e'^e tt ^ f • • •, where f — iorb quark in the form suitable for an application of 
the optimal observables (00). The most general (CP- violating anc? CP-conserving) 
couplings for jtt, Ztt and Wtb have been assumed. All fermion masses except rrit 
have been neglected and we have kept only terms linear in anomalous couplings. 
We have assumed the tagging efficiency at the level of 60% both for lepton and b 



for Pe- 
for Pe- 
for P,- 



0.2 and P^i 
0.2 and P^i 
0.2 and P,4 



0.2 
0.1 
0.1 



(4.15) 



for Pe- = 0.2 and Pe+ = 0.1 



quark detection, the range of the polar angle restricted by | cos^j.| < 0.9 and the 
integrated luminosity L — 500 fb~^. 

CP- violating charge forward-backward asymmetry A{;p has been introduced as 
an efficient way for testing CP- violation in top-quark couplings. Since the angular 
distribution for leptons is insensitive to variations of the standard V— A structure of 
the Wth coupling, the asymmetry could be utilized for a pure test of CP- violation 
in the top-quark production process. The expected statistical significance Ngr, for 
the measurement of the asymmetry has been calculated. We have found that it 
should be possible to detect A'ap ^-^^ (4.3(t) level for bottom quarks (leptons) for 
unpolarized beams, assuming CP- violating couplings of the order of 0.05. Having 
both beams polarized at 80% the signal for bottom quarks (leptons) could reach 
even 16(T (6.6cr). 

Next, the 00 procedure has been applied to the angular distributions. In 
the case of the lepton angular distribution, the expected statistical significance 
for signals of non-standard physics varies between 1.3(7 and 35.5(7 assuming non- 
standard form factors of the order of 0.05. It turned out that in the case of the 
bottom-quark angular distribution the statistical significance of the signal is in 
general higher than for leptons because of larger event rate and varies between 
1.8(7 and 85. lex for the same non-standard form factors. 

When deriving the above results we have fixed all the non-SM parameters to be 
-1-0.05 as a reasonable example for the strength of beyond-the-SM physics. How- 
ever, final results for statistical significances considered here depend on the size of 
the non-standard parameters. The most convenient beam polarizations for a mea- 
surement of the asymmetry A^^p and for testing the angular distributions varies 
with the non-standard parameters, as well. Therefore one should stress that the 
beam polarizations should be carefully adjusted for each model to be tested in 
actual experimental analysis. However, in any case, the above results show that a 
measurement of A^p and GO analysis of the angular distributions are both very 
efficient for new-physics search. 

Then we have analyzed the angular and energy distribution of the lepton toward 



separate determinations of the anomalous form factors. In order to reach the 
highest precision we have been adjusting beam polarizations to minimize errors 
for each form factor. We have found that at ^/s = 500 GeV with the integrated 
luminosity L = 500 fb^^ the best determined coupling would be the axial coupling 
of the Z boson with the error A[ Re(5Az) ] = 0.07 while the lowest precision is 
expected for Re((5D^) with Z\[Re((5D^) ] = 14.4. This result is independent of the 
choice of the non-SM parameters in contrast to the above two types of analyses. 

Concluding, we have observed that the angular distributions and the angular 
and energy distributions of top-quark decay products both provide very efficient 
tools for studying top-quark couplings to gauge bosons at linear colliders. 
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Appendix 



Integrals of Oi (x) denoted in the main text by i7/ in the angular distribution 



eq. 



are the following: 



= Z}W - (1 



+ 



f2i^*^ = 4 



2/32) - 2 Re(G'S*^) 
af[ 2(1 - p^)Re{D^y}) - Re{F}*^) + (3 - 2(3^)Re{G^^'>) 



D^y^ + ^i*^ + 2 

+ a^Re(2DW 
Re(4J) + 2«^[(l 



Re(Gi*^) 



,^--Re(FW-GW)] 



(*) 



1-/32 1 + P 



[]^^*^ = (3 - 2/32)[ + Z)W + 2Re(Gi*^) ] 



(*) 



+ 3a^[2(l-/32)Re(DW 



y^j-Re{Ft^) + {3-20')Re{G^;>)] 
+Z}J) +2Re(GS*^) 

+ «^Re(2DS - Fi*^ + 3Gl 



In 



1-/3' 



"^3 > 



(A.l) 



2/3 1-/3 

Next we present explicit formulas of the coefficient functions for the nine anoma- 
lous form factors in eq.([4.5|) J^^^^g cd}{7 ^) ^^"^ ^mi.'^^ ^) (/ = ^/^)- 

^(3 - /32)C(Dv : A,)f'{x) + 2a(,C{DvA : ^.)^?' (x) ] (1 + cos^ 
a^(l-/32)C(DvA:A)f (x) 

-^{C{Dv : A) + 2a^C(Z}yA : A,)}{2hi{x) -hi{x)}]{l-3 cos' 
+2\aiC{Ev:A,){gf{x) - hi{x)} + 2C{Eva: A,){ff{x) - hi{x)} ] cosO, 



(A.2) 



= ]^(3^C{Da : B,)f{x){3 - cos' 9) + 2a(,C{DvA : B,)gf{x){l + cos' 9) 

-i [ {C{Da : fi.) + 2ao'(l - [3^)C{Dva : A,)}/^ (x) 

-{C{Da : 5^) + 2a(^C{DvA ■ B^)}{2h{{x) - /i^(x)} ] (1 - 3 cos' 0; 
+2 [ {a^(l - /3')C(E4 : B,) + 2C{Eva : B,)}ff{x) + a',C{EA : 5,)(7^(x) 



-{aiC{EA : B,) + 2C(E,.^ : B,)}h{{x) 



cos 6*, 



(A.3) 



= -/3'C(Gi:a)/^(a;)(l + cos'^) 
+2a{iC{G2:C^)[f^'{x) + gf{x) - /i((x) ] cos^ 

{C(Gi : G,) + a5(2 - /3')C(G3 : C,)}/^ (x) + ^^^(Gs : a)(7^(x) 
-{2GiGi : a) + 3a^C(G3 : G,)}hi{x) 



+{CiGi:C,) + al,CiG3:C,)}hiix) (l-Scos^ 



(A.4) 



-r-/{*), 
Dv ' 



x,9) 



a(^C{Fi:D^)\ f^{x) - /i((a;) 1 (1 - 3 cos^ 9) - a(^C{Fi: D^)g^' {x){l + cos'^ 9) 



-2aiC{F^:D.,) f ^ {x) + {x) - h{{x) cos9, 



(A.5) 



while J-'2r\x, 9) takes different forms for f = i and f = b as 
1 r 



r,^n'{x,9) 



(3 - - (1 - 3p'')D'^'*> - 2(1 - 0')Re{D 



(0,*) 



+2Re(Dg*^)5(7^(x)(l + cos^ 



+ 2 

+2 



+ D)"'*^ + 2Re(DVT') 25/i^(x) - 5/i^(x) (1-3 cos' 



i(o,*) 



l(0,*) 



(1 - /3')4°'*) + 2Re(E|.^ 
) 5g\x) cos 9 



n(0,*)> 



+2 ( 4°'*) + 4°'*) 



-2 

1 r 
+ 2 



(0,*) 



+ E, 



(0,*) 



(3-/3')(Z}{.°'*)+DA'^'*0 + 6(l-/3')Re(Di.r) r(a:)cos'^, (A.6) 



2Re(E- 
i(o,*) 



(0,*)^ 



6f{x) COS 9 



Sh\{x) cos 9 



(0,*)^ 



and 



^5W(x,^) 

= {Re(Dgi*^) [ -{(1 - f3^)f{x) - 2h\{x) + /i2(x)}(l - 3 cos' ^) 

+ 2^''(x)(l + cos'0)" 
+2[ (1 - p')Ef*^nx) + (4°'*) + Ef*^){g\x) - h\ix)} ] cos9} , (A.7) 

where the functions f^{x), g^{x), h{ 2{x), Sf^{x), Sg^{x) and Sh{ 2{x) are defined as 

Ff{x) = r{x) + Re{f 2^)6 fix), 
Gf{x)=g^{x) + Re{f2^)6g^{x), 

HUx) = hl2{x) + Re(/2«)5M,2(x), (A.8) 
with F^{x), G^{x) and H(2{x) being given as follows [jl3| 



Ft dxduj ' 



G^{x) = — duo 
Bf J 



X- 



l+(3 



l-uo 



Fi dxdoj ' 



J duj{l — uj) 



1 (fTf 



Bf X J " ^ Ft dxdu ' 



X 



Ft dxdu 



(A.9) 



and oj is defined as a; = {pt — P/)^/""^? • 

After performing tlie above integrations using 



f 1 + 13 



Ft dxduj 



P W 



■UJ 



l + 2Re(/2^)x/F(^- 



uj 1 + 2r 



S{uj — r) 



for / = £+, 



for / = 6. 



I 2/5(1 -r) 

one obtains the following explicit forms of f^{x), g'{x), h{2{x), Sf'{x), Sg^{x) and 
Shi 2{x) for leptonic and bottom-quark final states: 
• For / = ^ 

3(1+^ 2 / 3(1+/?) , 2 



^^(x) = fix) + 



3(1 + py 



(3W 



X 



2I3W 
a; + In |1 — a;| 



h{{x) = 

K{x) = 



1-/3M 

2PW X 
1 



6^2(3 - 2uj) 



A(5W 



(l + /3)(l-/3)^- u;^(6-8a; + 3u;^) 



x^ 



5f\x) = - ^^Lt/^ V^[2a; + 21n|l-a;| 



5g\x) = r (^) 



6(1 + /?)^ 



l + 2r 



Ul- 



(3W 



X 



In |1 — a;| 



+ 



+ 



\-u 1 + 2r 



(u; + ln|l 



bh\{x) 



3(1 - 



/3iy X 



a;' 



,/ 3 — 2a;\ 
^ ~ l + 2r>' 



1 



2PW 



(1 + /?)(1-/?)^ 



3cc;2 



2a;^(3 - 2a;) - —(6 - 8a; + 3t^ 



(A.IO) 



where uj± are given as follows: 
For r > B (r = M^/m^ and B 



For r < B 



1 — r, 
1 — r, 
1 — X, 



For / 



r, 



1 





x/5 



= l-x/B 
= l-x/B 
= 

1 + , 



2/3(1 -r) ' 
1+P 



1 



-X 



for i?r < X < B 

for B < X < r 

for r < X < 1 

for i?r < X < r 

for r < X < B 

for i? < X < 1 

constant), 



(A.ll) 



(A.12) 



1-r ^2/?(l-r)^ 



2/?x ' 
(l-r)(l + /j)(l-/?)^ 

2/5x2 
: 5^'(x) = Sh\{x) 



6h\{x) = 0, 



(A.13) 



where x is bounded as 

B{1 — r) < X < 1 — r. 

The coefficients C {X : Y) employed in the definition of the coefficient functions 
have been introduced through the following formulas: 



v=-y,Z 



i(o,*) 



Re(Dg: 



(0,*)^ 



Re(DVTO 

+ 5] [C(Z}yA:^.)Re(M,) + C(Dy^:5,)Re(5B,)], (A.14) 

D=7,Z 



and in the analogous manner for Ey^AyA-, -^1-4 and (71^4. Dy'^yj^^ _E'y^''^, Fil,4 , 



i(0,*) 



.(0,*) 



G^I:Ca could be obtained from eq. (|A.17|) below as a SM approximation of D'y' 



V,AyAi 



^v^a,vaj -^1^4) ^^1^4- Explicit forms of the independent coefficients are given as 

C{Dv : A) = 2C[V^A^ - (Pe + VeV^)d'Az ] , 
C{Ev:A^) = -2C[ V^A^ - {V^ + v,V^)d'Az ], 
C{Dva:A^) = -C{Ve + VeV^)d'Bz, 

C{Eva:A^) = C{V^ + v,V^)d'Bz, 

C{Dv:Az) = -2C[ {Ve + v,V^)d'A^ - {2v,V^ + (1 + vl)V^}d'^Az ], 
C{Ev:Az) = 2C[ (V^ + VeVe)d'A^ - {Iv^V^ + (1 + vl)V<^}d'''Az ], 

C{Dva:Az) = C[ 2v,V^ + (1 + vl)V^ ]d'^Bz, 

C{Eva:Az) = -C[ 2veV^ + (1 + vl)V^ ]d'^Bz, (A.15) 

where Vg = —l+Asin^OiY, d' = s/[4sin6'vi/ cos6'vi/(s— M|)], two polarization factors 
and Pig, are defined as 

= Pe- + P^ = 1 + A- A+, 

and the others are thereby given as 

C(Da:B,) = 2C(Dva:A,), C(Ea:B,) = 2C(Eva:A,), 
C{Dva:B,)^C{Dv:A,)/2, C{Eva:B,) ^ C{Ev:A,)/2, 



(A.16) 



C(Gi : = C{Dv ■.A,)/2, C{G2 : C,) - C(£;y : A)/2, 

C(G'3:a) = ^(6-4: a) = C{Eva:A,), 

C{F, : D,) = -CiDv : A)/2, C(P2 : D,) = -C{Ev : A„)/2, 

C{Fs:D,) = -C{Dva:A,), C{F^:D,) = -C{Eva:A,). 

As explained in the main text, they are not always independent of each other. 
When Pq = iP®, i.e., Pg- = Pe+ = ±1, we have 

C{{Dv, Ev, DvA, Eva}:Az) = t(1 ± v,)d'C{{Dv, Ey, Dva, Eva}:A^), 

As a consequence of the above relations one gets 



In this case all we can determine (for the production form factors) are the following 
four combinations 



Re{6{A, B, C, D}^ ^ (1 ± Ve)d'6{A, B, C, D}z). 

Finally we present here formulas for -D[^^_^, Ey^j^^yj^^, -fi~4, GiIa foi' complete- 
ness: 



r)(*) 

^V,A,VA 


= 'P^Dv,A,VA 


— '^e-E'y, A, VA 


rni*) 

^V, A, VA 


— 'P^Ev,A,VA - 


- V®Dv,A,VA; 


-^1,2,3,4 ■ 


= 'P(g)-^1,2,3,4 - 


■ "^0-^2, 1,4, 3, 


<-^l,2,3,4 


= 7^0(^1,2,3,4 - 


- P©G2,1,4,3, 



-2{A^Ved' - Az{\ + t;f)rf''}Re(Mz) ], 
Da = C[Bl{l + vl)d'^-2BzVed'Re{5B^) + 2Bz{l + vl)d'^Re{5Bz)i 
DvA = C[ -A^BzVed' + AzBz{l + vl)d'^ - BzVed'{SA^)* 

+ {A^ - Ved'Az)SB^ + Bz{l + vl)d'\5Azy 

-{A^ved' - Az{l + vl)d'^}6Bz], 
Ev = 2C[A^Azd' -Alved'^ + Azd'Re(5A^) + (A^d' -2AzVed'^)Re(5Az)], 
Ea = 2C[ -Blved'^ + Bzd'Re{6B^) - 2BzVed'^Re{6Bz) ], 
Eva = C [A^Bzd' - 2AzBzVed''^ + Bzd'{6A^y + Azd'6B^ 

-2BzVed'\5AzT + {A^d' - 2AzvJ'^)5Bz], 
Fi = C[-{A^-AzVed')6D^ + {A^Ved'-Az{l + vl)d'^}SDz], 
F2 = C[ -Azd'SD^ - {A^d' - 2AzVed''^)6Dz], 
Fs = C[ BzVed'SD^ - Bz{l + vl)d'^SDz], 
F4 = C[-Bzd'5D^ + 2BzVed'^5Dz], 

Gi = C[{A^-AzVed')SC^-{A^Ved'-Az{l+vl)d'^}6Cz], 



G2 = C[Azd'6Cy + {A^d' - 2AzVed'^)5Czl 
G3 = C[-BzVed'5C^ + Bz{l+vl)d'HCzl 

Gi = C[ Bzd'SG^ - 2BzVed'^6Cz ] (A. 18) 

with C = 1/(4 sin^ Ow)- 
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